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Abstract. We detect near-infrared H2 emission in DG Tau using the Infrared Camera and Spectrograph (IRCS) 
on the 8.2-m SUBARU telescope. The spectra obtained along the jet axis show that the centroidal position of the 
1-0 S(l) emission is offset by 0.2" from the star towards the jet, while those obtained perpendicular to the jet axis 
show a marginal extension, indicating that the emission line region has a typical width of ~0.6" . Their line profiles 
^v^j show a peak velocity of ~15 km s _1 blueshifted from the systemic velocity. These results indicate that the emission 

originates from a warm molecular wind with a flow length and width of ~40 and ~80 AU, respectively. The line 
flux ratios (Ii-os(o)/Ii-os(i) an d an upper limit for /2-is(i)/^i-os(i)) suggest that the flow is thermalized at a 
temperature of ~2000 K, and is likely heated by shocks or ambipolar diffusion. The observed velocity and spatial 
extension suggest that the H2 and forbidden line emission originate from different components of the same flow, 
i.e., a fast and partially ionised component near the axis and a slow molecular component surrounding it. Such a 
pH I flow geometry agrees with model predictions of magneto-centrifugal driven winds. 
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' ° ular envelope (Wiseman et al. 2001; Testi et al. 2002) 

Among the models proposed to date, magneto 
centrifugal wind models seem the most promising as the 
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limited fraction of the entire flow close to the jet axis (e.g., 
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ponent has been suggested by the presence of a cavity in 
the molecular core/envelope (see Davis ct al. 2002, and 
references therein), and a broad range of morphology of 
molecular bipolar outflows (see e.g., Lee et al. 2000, 2001). 
In addition, Davis et al. (2002) and Saucedo et al. (2003) 
revealed the presence of shock-heated H2 along the cavity 
wall in some young stellar objects (YSOs), presumably due 
to interaction with a wide-angled wind. Direct detection 
of such a flow component is desired to test the validity of 
the magneto-centrifugal wind scenario, to investigate the 
angular transfer problem in detail, and to understand the 
nature of molecular outflows. 

We report the detection of near-infrared H2 emission 
in DG Tau, which originates from a warm molecular wind 
close to the star. The object is one of the most active T 
Tauri stars known, and exhibits an energetic jet like Class 
I protostars (e.g., Eisloffel & Mundt 1998; Davis et al. 
2003). Since the object is not heavily embedded like Class 
I protostars (Av=2.2 — Muzerolle, Calvet, & Hartmann 
1998), it has been observed as one of the best targets for 
studying outflowing activities close to the star (see, e.g., 
Solf & Bohm 1993; Kepner et al. 1993; Bacciotti et al. 
2000, 2002; Pyo ct al. 2003). In particular, Bacciotti ct al. 
(2000, 2002) observed the forbidden line outflow using the 
Hubble Space Telescope, and revealed the presence of an 
"onion-like" kinematic structure, i.e., a continuous brack- 
eting of the higher speed central flow within the lower 
speed, less collimated, broader flow. 

In this paper, we show that the warm molecular wind 
appears as a part of the "extensive" onion-like kinematic 
structure observed in the forbidden line outflow. The rest 
of the paper is organised as follows: we describe the de- 
tails of the observations and results in Sect. 2; we compare 
our results with previous observations of near-infrared H2 
emission in other YSOs, and investigate the driving, lo- 
cation of the launching region, and excitation/heating of 
the warm molecular wind in Sect. 3; we then give our 
conclusions in Sect. 4. Throughout the paper, we adopt a 
distance to the object of 140 pc based on HIPPARCOS 
observations (Wichmann et al. 1998). 

2. Observations and Results 

Observations were made on the night of 2002 November 
25 at the SUBARU 8.2-m telescope using the Infrared 
Camera and Spectrograph IRCS (Tokunaga et al. 1998; 
Kobayashi et al. 2000). The echelle grating mode with 
a 0.3" wide slit provides a spectral resolution of ~10 4 . 
The pixel scale of 0.06" provides good sampling of the 
seeing profile (0.7-0.8" during the observations), thereby 
allowing us to investigate the spatial structure of the 
emission line region at sub-arcsecond scales. The spectra 
were obtained along, and perpendicular to, the jet axis 
(P.A.=226° — Bacciotti et al. 2000) with a total integra- 
tion of 480 s for each direction. The cross disperser was 
set to observe the following wavelength coverages simulta- 
neously: 1.958-2.011 (ordcr=28), 2.031-2.085 (27), 2.109- 
2.165 (26), 2.193-2.251 (25), 2.285-2.344 (24), and 2.384- 



2.446 /im (23). In addition to the target, A-type bright 
standards were observed at similar airmasses to correct 
for telluric absorption. The flat fields were made by com- 
bining many exposures of the spectrograph illuminated by 
a halogen lamp. 

The data were reduced using the KAPPA and 
FIGARO packages provided by Starlink. We obtained 
position-velocity diagrams via standard reduction pro- 
cesses: i.e., dark-subtraction, flat-fielding, removal of bad 
pixels, correcting for curvature in the echelle spectra, 
wavelength calibration, correcting for telluric absorption, 
and night-sky subtraction. Wavelength calibrations were 
made by identifying OH airglow and telluric absorption 
features. For the templates, we used the UKIRT OH line 
list 1 and model spectra provided by ATRAN (Lord 1992). 
The accuracy of wavelength calibration is better than half 
a pixel on the detector, corresponding to 4 km s _1 . The 
systemic motion of the object was calibrated based on pre- 
vious observations of Li I A6707 photospheric absorption 
(U Hc i=16.5 km s^ 1 , Bacciotti ct al. 2000). 

We detected H 2 1-0 S(l) (2.1218 fim), 1-0 S(0) (2.2233 
/im), 1-0 Q(l) (2.4066 ^m) and 1-0 Q(3) (2.4237 fim) 
emission in our cross-dispersed spectra. Since the latter 
two lines severely suffer from telluric absorption, we use 
the former two lines and the upper limit of the 2-1 S(l) 
(2.2477 fim) flux for the rest of the paper. Fig. 1 shows 
the position- velocity (P-V) diagram of the 1-0 S(l) emis- 
sion, and its spatial distribution integrated over —40 to 
10 km s _1 . Although we did not fully resolve the emission 
line region, 1-0 S(l) emission still exhibits a distribution 
different from the continuum emission. In the spectra ob- 
tained along the jet, the position of the 1-0 S(l) emission is 
displaced from the continuum position by ~0.2" towards 
the jet. Those perpendicular to the jet show broader and 
slightly asymmetric distribution about the continuum. To 
investigate the spatial scale of the emission line region 
here, we convolved the seeing profile at the continuum 
with a Gaussian, and plot it in Fig. 2 together with the 
distribution of the 1-0 S(l) emission. The figure shows that 
1-0 S(l) emission has a typical angular scale of ~0.6" per- 
pendicular to the jet. The 1-0 S(0) line is much fainter than 

1- S(l), and its spatial distribution is the same as that 
of 1-0 S(l) within the uncertainty of the measurement. 

Fig. 3 shows the spectra at 1-0 S(l), 1-0 S(0), and 

2- 1 S(l), integrated over ±1.2" along the jet and cen- 
tred on the star. The lines exhibit a peak at — 15±4 km 
s _1 from the systemic velocity. We believe this is not 
a spurious wavelength shift due to uneven illumination 
of the slit (cf. Bacciotti et al. 2002), since (1) the se- 
lected slit width of 0.3" is much smaller than the see- 
ing (0.7"-0.8"), and (2) the spectra at different slit posi- 
tions show the same velocity within the uncertainty of the 
wavelength calibration. The FWHM of 35 km s _1 in the 
1-0 S(l) profile is comparable to the spectral resolution 
(~30 km s _1 ), implying that the internal velocity disper- 

1 http://www.jach.hawaii.edu/JACpublic/UKIRT/astronomy 
/calib/oh.html 
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Fig. 1. (top) Continuum-subtracted position-velocity diagrams of 1-0 S(l) emission along and perpendicular to the 
jet axis. The contour spacings are 10% of the peak intensity. The subtracted continuum flux is measured closest to 
the dot-dashed lines. The continuum flux is not constant along the wavelength due to a number of shallow absorption 
features, and this causes residual after continuum subtraction at ^100 km s _1 . (bottom) Intensity distribution of 1-0 
S(l) integrated between dot-dashed lines in the P-V diagrams. The grey histogram shows the seeing profile observed 
at the continuum. 
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Fig. 2. Same as the bottom right of Fig. 1, but the seeing 
profile measured at the continuum is convolved with a 
Gaussian having a FWHM of 0.6" (grey solid), 0.3" and 
0.9" (grey dashed), respectively. The centroidal position 
of the convolved profiles are adjusted to fit the 1-0 S(l) 
intensity emission (black solid). 



sion is much smaller. Assuming a K magnitude of 6.74 
(Muzerolle et al. 1998), we derive the 1-0 S(l) and 1-0 
S(0) fluxes of 1.8±0.1xl0" 14 and 3.7±0.5x 10~ 15 erg s" 1 
cm~ 2 , respectively, and a 3-er upper limit of the 2-1 S(l) 
flux of 2.2xl0~ 15 erg s" 1 cm" 2 . The estimated 1-0 S(l) 



flux is consistent with the upper limit for the 1-0 S(l) flux 
(<5.3 x 10~ 14 erg s" 1 cm" 2 ) obtained by Carr (1990). 

Table 1 shows the 1-0 S(0)/l-0 S(l) and 2-1 S(l)/1- 
S(l) ratios measured in the same region as for Fig. 3, 
i.e., within ±1.2" along the jet axis centred on the star. 
For comparison, we also tabulate the line ratios for lo- 
cal thermal equilibrium (LTE) at T =1000-3000 K, where 
we have calculated these using the ^-coefficients provided 
by Turner, Kirby-Docken, & Dalgarno (1977). The table 
shows that observed line ratios are consistent with ther- 
mal excitation with a temperature of ^2000 K. Such a 
spectrum is often observed in shocks in extended outflows 
and Herbig-Haro objects (see e.g., Gredel 1994; Eisloffel, 
Smith, & Davis 2000b), although this could also result 
from thermal collisions in a dense UV-photodissociation 
region (see e.g., Sternberg & Dalgarno 1989; Burton, 
Hollenbach, & Tielens 1990) or via X-ray heating (see e.g., 
Maloney, Hollenbach, & Tielens 1996; Tine et al. 1997). 
We investigate this in detail in Sect. 3.3. 

The blueshifted line profile and spatial offset away 
from the star indicate that the detected H 2 emission orig- 
inates from the outflowing gas. The absence of red-shifted 
components is attributed to the obscuration of the coun- 
tcrflow by a circumstellar disk, common for forbidden line 
emission (e.g., Appenzeller, Jankovics, & Ostreicher 1984; 
Edwards et al. 1987). The radial velocity of —15 km s _1 
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Table 2. H 2 Outflow Parameters 
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Fig. 3. Spectra at 1-0 S(l), 1-0 S(0), and 2-1 S(l) ob- 
tained along the jet. The P-V diagram obtained along the 
jet is integrated over ±1.2" of the continuum position, and 
normalized to the continuum at 2.12 firm. The dot-dashed 
lines indicate the velocity range in which we measure the 
line flux. The asterisks indicate photospheric absorption 
features. Relative line fluxes between the spectra are con- 
served since the spectrograph covers all the lines with a 
single exposure. 

Table 1. H 2 Line Ratios" 

observed calculated 

(T=1000 2000 3000 K) 



h 
h 



-is(i)/h- 
-os(0)/h- 



0S(1) 
03(1) 



<0.11° 
0.20±0.03 



0.005 
0.269 



0.082 
0.211 



0.210 
0.194 



Extinction is not corrected. 
6 An upper limit for the 3-<r level. 



is slower than the forbidden line outflow, which exhibits 
two blucshifted peaks at —200 and —20-100 km s _1 (e.g., 
Edwards et al. 1987; Hamann 1994; Hartigan et al. 1995; 
Takami et al. 2002; Pyo et al. 2003). This trend is also ob- 
served in some Class I protostars, as shown by Pyo et al. 
(2002) and Davis et al. (2003). The measured flow width of 
—0.6" in H2 emission is wider than the forbidden line out- 
flow, which has FWHM widths of <0.1" and -0.2" in high 
and low velocity components, respectively, at the same dis- 
tance from the star (Bacciotti et al. 2002). These results 
indicate that H 2 emission originates from a different (i.e., 
slower and wider) flow component to the forbidden line 
outflow. 

Table 2 summarises the flow parameters of the warm 
H 2 wind. The wind is not fully resolved in our spectra, 
thus the flow scale along the jet axis is estimated using 
the measured spatial offset of 0.2" , and assuming an incli- 
nation angle the same as the jet (45° — Pyo et al. 2003). 
The derived spatial scale of —40 AU is about half as large 
as that perpendicular to the jet axis (—80 AU), indicat- 
ing that the flow has an opening angle of -90° if we as- 
sume zero width at the base. The luminosity is obtained 
assuming the extinction of Ak=0.2, based on the visual 
extinction to the object of Ay=2.2 (Muzerolle et al. 1998) 



Mean velocity" 

Flow scale along the jet axis a 

Flow scale perpendicular to the jet axis 

1-0 S(l) luminosity 

H 2 mass emitting 1-0 S(l) 

Total H2 mass 6 

Mean H2 number density 6 

Mass loss rate 6 

Momentum flux 6 

Kinematic luminosity 6 (M H2 V :i /2R) 



-20 km s _i 
-40 AU 
-80 AU 
1.3 x 10~ 5 L Q 
2.7 x 10~ 10 M Q 

2.1 x 10~ 8 M© 
^ 4xl0 4 cm~ 3 
^ 2.2 x 10~ 9 M yr" 1 
^ 4.4 x 10~ s M km s" 
J> 7.2 x 10~ 5 L© 



a Corrected for inclination to the line of sight. 
6 The lower limit is derived assuming the local thermal 



equilibrium at T =2000 K. See text for details. 



and A K ^0.1A V (Becklin et al. 1978). The total H 2 mass 
is estimated assuming Mi_ s(i)/Mh 2 = 1.28 x 10~ 2 at 
T = 2000 K, based on Smith & Mac Low (1997) and 
Davis et al. (2001 ) 2 . This mass should be a lower limit 
since the flow may include colder gas, which does not con- 
tribute to ro- vibrational H 2 emission. Indeed, the presence 
of such cold gas is predicted by shock (see Eisloffel et al. 
2000b, and references therein) and UV-photodissociation 
models (e.g., Burton et al. 1990). The mean H 2 number 
density is estimated assuming a conical geometry with a 
filling factor of 1. This value should also be a lower limit 
due to the underestimated H 2 mass, and the fact that the 
filling factor may be less than 1 due to a clumpy or hollow 
geometry. 



3. Discussion 

3.1. Comparison with H 2 emission in other YSOs 

Unlike CO and CS, the H 2 molecule has no net dipole 
moment, hence it is not readily observed in ambient cold 
(<ICl00 K) molecular clouds. However, shocks or UV/X- 
ray radiation in YSO allows the molecule to be excited to 
ro-vibrational states, thereby allowing us to probe active 
circumstellar regions. 

In the last decades, a number of molecular outflows 
and Hcrbig-Haro objects have been found to exhibit near- 
infrared H 2 emission (see., e.g., Gautier et al. 1976; Grcdcl 
1994; Davis et al. 1998; Chrysostomou et al 2000; Eisloffel 
et a. 2000b). Their morphology, line profiles, and line 
flux ratios are well explained by thermal excitation from 
shocks. In addition to the extended jets/outflows, their 
driving sources (Class I protostars) also exhibit near- 
infrared H 2 emission (Carr 1990; Greene & Lada 1996; 
Reipurth & Aspin 1997). Recent high-resolution spec- 
troscopy and Fabry-Perot observations have confirmed 
that the emission originates from the base of the jet 



2 Assuming a population fraction Xv,j(T)=g v ,jZ~ 1 (T)x 
exp(—AE v ,j/kT), where g Vt j is the statistical weight, AE V; j 
is the energy difference from the ground state, Z(T) is the 
partition function Z(T)=0.024T/(l-exp(-6000/T)), and T is 
the temperature. 
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(Molecular Hydrogen Emission-Line Regions — Davis et 
al. 2001, 2002). Their line profiles show properties similar 
to forbidden line emission in T Tauri stars: (1) the profiles 
often contain low (5-20 km s _1 ) and high velocity (50- 
150 km s _1 ) components, (2) the low- velocity component 
(LVC) is more common than the high- velocity component 
(HVC), and (3) the HVC is further offset than the LVC. 

In addition to Class I protostars, several T Tauri stars 
are known to exhibit ro-vibrational H2 emission. T Tau 
exhibits bright H 2 emission up to ^1000 AU scale, and 
spectro-imaging by Herbst et al. (1996, 1997) show that 
the emission originates through interaction between out- 
flowing and ambient gas. In contrast, H 2 emission in other 
T Tauri stars is likely to be associated with quiescent 
gas in the circumstellar disk. Weintraub, Kastner, & Bary 
(2000) and Bary, Weintraub, & Kastner (2002, 2003) show 
that the 1-0 S(l) line profile in these objects are cen- 
tered at the star's systemic velocity, in contrast to the 
blucshifted emission of Class I protostars. The observed 
line widths of 9-14 km s _1 and a double-peak profile in one 
of the objects (LkCa 15) are well explained by Keplerian 
motion in the circumstellar disk. The authors suggest that 
the H 2 molecules in these objects are excited by X-ray or 
UV radiation. 

As shown in Sect. 2, near- infrared H 2 emission in DG 
Tau is associated with outflowing gas within 100 AU of 
the star, sharing the same origin as Class I protostars 
rather than T Tauri stars. This result is attributed to 
its evolutionary phase. Although DG Tau exhibits opti- 
cal properties similar to T Tauri stars, its flat spectral 
energy distribution (Adams, Emerson, & Fuller 1990) and 
energetic jet (e.g., Eisloffel & Mundt 1998; Dougados et 
al. 2000; Davis et al. 2003) indicate that the object is in a 
transient phase between a Class I protostar and a T Tauri 
star (see, e.g., Greene & Lada 1996 for discussion). The 
1-0 S(l) luminosity in the DG Tau wind is 10-100 times 
fainter than Class I protostars (typically 10 -4 -10~ 3 L Q 
— Davis et al. 2002), and unlike some Class I protostars, 
a high- velocity component is absent. This could simply be 
a consequence of the evolution of DG Tau from the Class 
I to the Class II phase. 

Bary et al. (2003) detected quiescent 1-0 S(l) emission 
in GG Tau A, TW Hya, LkCa 15, and DoAr 21, and the 
measured fluxes in these objects correspond to luminosi- 
ties of 4.2xl0~ 6 , I.lxl0- 7 ,'l.0xl0- 6 , and 9.1xl0~ 6 L Q , 
adopting distances of 140, 60, 140, and 140 pc, respectively 
(Wichmann et al. 1998; de Geus & Burton 1991). These 
luminosities are smaller than the measured 1-0 S(l) lumi- 
nosity in the DG Tau wind, typically at least by a factor 
of 3. Thus, 1-0 S(l) emission from the quiescent gas could 
contribute to the total 1-0 S(l) luminosity in DG Tau, but 
the total 1-0 S(l) luminosity should still be dominated by 
emission from the wind, consistent with the observed line 
profiles and offset shown in Sect. 2. 



3.2. Driving mechanism and wind launching region 

Observations of a wide-angled wind will allow us to inves- 
tigate the driving mechanism of the jet/ wind in detail. In 
particular, the low extinction of DG Tau allows us to com- 
pare the kinematics of the wind with the optical forbidden 
outflow, thereby allowing us to investigate the entire flow 
structure close to the driving source. 

The Hubble Space Telescope observations by Bacciotti 
et al. (2000, 2002) revealed the detailed kinematics in the 
DG Tau jet with a resolution of ~0.1" . The authors show 
that the centroidal velocity of the low velocity compo- 
nent monotonically decreases with distance from the jet 
axis, exhibiting an "onion-like" kinematic structure. At 
0.2" from the star, the component shows a centroidal ve- 
locity of ~60±20 km s _1 at the axis, and ^30±10 km s 
0.2" from the axis (sec Fig. 1 of Bacciotti et al. 2002). 
If we extrapolate the velocity field in the forbidden line 
jet to the typical width of the warm H 2 wind (i.e., ±0.3" 
from the axis), we could expect a velocity of ~15 km s _1 , 
in excellent agreement with the measured velocity in H 2 
emission. This indicates that the warm molecular wind in 
H 2 emission can be attributed to an outer extension of 
the onion-like kinematic structure observed in the forbid- 
den line outflow. This agrees with the magneto-centrifugal 
wind models, which predict that (1) fast collimated com- 
ponents and slow poorly-collimated components coexists 
(e.g., Shu et al. 2000), and (2) forbidden line emission 
traces only a limited fraction of the entire flow close to 
the jet axis (Li & Shu 1996; Cabrit et al. 1999; Shang et 
al. 2002). 

In addition to the magneto-centrifugal wind scenario, 
the presence of the warm molecular wind may also be ex- 
plained by entrainment of the surrounding material by 
the collimated jet. This mechanism has been proposed, 
for instance, for H 2 emission in Class I protostars (Davis 
et al. 2003) and the low velocity component of forbid- 
den line emission (Pyo et al. 2003). Although our data do 
not exclude this possibility, it may have a few problems. 
Model calculations by Li & Shu (1996) suggest that such 
interaction close to the driving source is weak because of 
"magnetic cushioning", i.e., the effects of frozen-in mag- 
netic field in the flow. In addition, 13 CO observations by 
Kitamura, Kawabe, & Saito (1996) suggest that the ma- 
jority of the matter surrounding the jet has already been 
blown away, leaving nothing to entrain. These authors also 
show that the remnant, flat envelope at a ~1000 AU scale 
is slowly (~1.5 km s _1 ) expanding, and it may be diffi- 
cult to understand such motion without invoking with a 
wide-angled wind. 

Since the majority of the surrounding matter has al- 
ready been blown away (Kitamura et al. 1996), the molec- 
ular wind must emerge from the circumstellar disk, i.e., 
the remaining reservoir of the H 2 molecules. This fact may 
give us a clue to the location of the jet/wind launching re- 
gion, the position of which have been debated for over 
a decade. Magneto-centrifugal wind models suggest two 
possibilities: (1) the X-wind model predicts that the flow 
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emerges between the inner disk edge and stellar magneto- 
sphere within a few stellar radii (Shu et al. 2000); and (2) 
the disk wind model predicts that the flow emerges from a 
wider (>1 AU) disk surface (e.g., Konigl & Pudritz 2000; 
Ferreira 1997). The disk wind model can easily explain 
the presence of H2 molecules in the wind, whereas it is 
not clear whether the X-wind model can do so. This is be- 
cause any H 2 molecules are likely to become dissociated 
close (<0.1 AU) to the star because of the harsh environ- 
ment. Further modelling including molecular chemistry is 
necessary to investigate this in detail. 

3.3. Excitation and heating 

Understanding the heating in the jet/wind is desirable, in 
particular in close proximity to the driving source, since it 
will allow us to perform comparisons with kinematic mod- 
els and observations, and investigate their driving mech- 
anism in detail. In the case of the collimated jet, several 
heating mechanisms have been proposed, including shocks 
(e.g., Hartigan et al. 1995), turbulent dissipation in a vis- 
cous mixing layer (e.g., Binette et al. 1999), ambipolar dif- 
fusion (Saner 1993; Garcia et al. 2001), and X-rays from 
the star (Shang et al. 2002). Recent high-resolution obser- 
vations (e.g., Bacciotti et al. 2000; Dougados et al. 2000) 
and detailed analysis of optical forbidden lines (Lavalley- 
Fouquet, Cabrit, & Dougados et al. 2000; Dougados, 
Cabrit, & Lavalley-Fouquet 2002) indicate that shocks in 
the form of internal working surfaces are likely to be the 
cause of heating on 10-100 AU scales. This conclusion is 
corroborated by the detection of Hcl emission in the DG 
Tau jet (Takami et al. 2002). On the other hand, Takami et 
al. (2001, 2003) revealed the presence of Ha jets on a few 
AU scales from some T Tauri stars. Based on their spa- 
tial scale and line profile, the presence of another heating 
mechanism in the inner (<10 AU) region is suggested. 

What heats and excites molecular hydrogen in the 
warm wind? The cooling timescale of molecular gas at 
~2000 K is less than a year (see Hollenbach & Natta 1995), 
smaller than the dynamical scale of the warm molecular 
wind we observed (~10 yr). This indicates the presence of 
a heating mechanism, which keeps the wind warm and al- 
lows H2 emission to be observed. Possible mechanisms in- 
clude: (1) shocks, which is responsible for near-infrared H 2 
emission in the extended jet/outflows (e.g., Gredel 1994; 
Davis et al. 1998; Chrysostomou et al 2000; Eisloffel et 
al. 2000b); (2) UV fluorescence, which is responsible for 
UV H 2 emission in accreting T Tauri stars (Valenti, Johns- 
Krull, & Linsky 2000; Ardila et al. 2002); (3) X-ray, which 
may be responsible for near-infrared H2 emission in the 
circumstellar disk of T Tauri stars (Weintraub et al. 2001; 
Bary et al. 2002, 2003); and (4) ambipolar diffusion (Safier 
1993). 

If the warm molecular wind consisted of a dense UV- 
photodissociation region (PDR), we could expect a 1-0 
S(l) surface brightness of up to ^10~ 3 erg s" 1 cm~ 2 
str -1 based on (1) the UV continuum flux from the 



star (~2xl0~ 13 erg s" 1 cm" 2 A" 1 at 1400-2000 A - 
Gullbring et al. 2000), (2) the spatial scale of the warm 
molecular wind, and (3) PDR models by Burton et al. 
(1990). This would require a surface area of at least 
— 4xl0 30 cm 2 to account for the observed 1-0 S(l) lumi- 
nosity. Together with a typical column density of the H 2 
emitting region (~T0 21 cm ~ 2 — Burton et al. 1990) and 
the flow velocity shown in Table 2, we derive a lower limit 
of the momentum flux of 7xl0 -6 M Q km s" 1 yr -1 . This 
flux is extremely large for the evolutionary phase of this 
object: it is even much larger than a typical momentum 
flux of molecular outflows associated with its younger pro- 
genitors, i.e., Class I protostars (2xl0~ 6 M Q km s _1 yr -1 
Bontemps et al. 1996). We then conclude that UV exci- 
tation is not likely for the heating of the warm molecular 
wind. 

X-ray photodissociation region (XDRs) models by 
Maloney et al. (1996) show that the X-ray surface bright- 
ness is a function of the X-ray luminosity at the radiation 
source, distance to the XDRs, column density, and hydro- 
gen nucleus density. Based on their models, we derive a 
1-0 S(l) surface brightness of 10~ 6 erg s _1 cm~ 2 str -1 
providing (1) a typical X-ray luminosity of a classical T 
Tauri star of 10 29 erg s" 1 (Neuhauser et al. 1995); (2) hy- 
drogen nucleus density of 10 5 cm~ 3 ; (3) a column density 
of 10 21 cm~ 2 , the same density as assumed for UV ex- 
citation described above. The derived surface brightness 
described above is smaller than that for UV excitation 
by a factor of ^10 3 , implying that a larger surface area 
or column density would be required to account for the 
observed 1-0 S(l) luminosity. However, such assumptions 
would provide even larger flow momentum flux than we es- 
timated for UV excitation. We conclude that neither UV 
or X-rays constitute the primary mechanisms. Note that 
the hydrogen nucleus density in the warm molecular wind 
may be higher than we assumed above (see Table 2), al- 
though a higher density would not increase the 1-0 S(l) 
luminosity. This is because it would cause a larger cooling 
rate, thereby allowing X-ray heating to be less effective. 

In contrast, shocks are a likely mechanism for the heat- 
ing of the warm molecular wind, since they are observed 
in the jet/outflows of various YSOs. The flow velocity of 
~20 km s _1 in the DG Tau wind is much larger than the 
sound velocity (3 km s _1 at T=2000 K), indicating that 
shocks could occur in the same manner as for extended 
jets and outflows: i.e., as internal working surfaces (e.g., 
Chrysostomou et al. 2000; Davis et al. 2000) or super- 
sonic turbulence (Davis & Smith 1996). In addition, the 
observed line flux ratios shown in Table 1 are similar to 
those in shock-excited extended outflows, which indicate 
thermal temperatures of 1000-4000 K (see., e.g., Gredel 
1994; Eisloffel et a. 2000b), supporting the shock heating 
scenario. 

In addition to shocks, ambipolar diffusion may also ex- 
plain the heating in the warm molecular wind we observe. 
Safier (1993) shows that H 2 molecules could survive in the 
outer region (>1 AU from the flow axis) of an onion- like 
kinematic structure, and they could be heated up to a few 
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thousands of Kelvins, consistent with our results. He also 
predicts that the temperature in the wind increases with 
the distance. The inner structure of the observed wind 
could be resolved by use of adaptive optics on 8-10 m 
telescopes, and such observations in the near future would 
allow us to identify the heating mechanisms, testing the 
validity of existing models. 

4. Conclusions 

We detected near-infrared H 2 emission in DG Tau, one of 
the most active T Tauri stars known. The spatial extension 
along and perpendicular to the jet, and measured blucshift 
of 15 km s _1 , indicate that the emission originates from 
a warm molecular wind with a flow length and width of 
^40 and ~80 AU, respectively. Via detailed comparison 
of the spatial scale and velocity field, we conclude that 
H2 and forbidden line emission originate from different 
regions of the entire flow, which includes (1) fast, well- 
collimated and partially ionised streamlines, and (2) slow, 
poorly-collimated, and molecular streamlines. Such a flow 
geometry agrees with predictions of magneto-centrifugal 
disk wind models. Although it could also be explained 
by the entrainment of the surrounding material by the 
collimated jet, we consider it unlikely. 

The measured line flux ratios (/i_os(o)/^i-os(i) an d 
the upper limit of l2-is(i)/Ii-os(i)) indicate that the flow 
is thermalized at a temperature of ^2000 K. The flow is 
likely to be heated by shocks or ambipolar diffusion. 
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